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We describe and analyse numerically schemes (i) for entan- 
gling orthogonal motional modes of one or a few harmonically- 
trapped atoms or ions, and (ii) for transferring the entan- 
glement from one of these local modes to a distant trapped 
atom (or atoms) via a light-mediated quantum state trans- 
fer procedure proposed in previous work [A.S. Parkins and 
H.J. Kimble, J. Opt. B: Quantum Semiclass. Opt. 1, 496 
(1999)]. Possibilities arising from these schemes include the 
generation of an Einstein-Podolsky-Rosen state in the posi- 
tions and momenta of distantly-separated trapped atoms and 
the preparation of delocalized mesoscopic vibrational states. 

PACS numbers: 03.67.Hk, 32.80.Lg, 42.50.-p 



I. INTRODUCTION 

In addition to being of fundamental interest in physics, 
quantum entanglement is an essential resource for quan- 
tum information processing and distribution. Popular 
candidates for experimental investigation in this con- 
text include trapped atoms, cavity quantum electrody- 
namics (cavity QED), and nonclassical light fields, using 
vifhich a variety of impressive quantum state manipula- 
tions have now been demonstrated, including elementary 
quantum logic operations 0-§Li preparation of quantum- 
mechanically entangled pairs 1 1- ^ and quadruplets j|] of 
atoms, and quantum teleportation |p|-|rT[|. 

In recent work we have proposed a scheme that offers 
the possibility of combining the features of trapped atoms 
and nonclassical light fields in a distributed network ; 
in particular, a scheme that enables motional quantum 
states to be coupled to propagating (nonclassical) light 
fields via interactions in cavity QED. This coupling en- 
ables the deterministic generation and distribution of 
quantum entanglement between different atoms and/or 
different light fields. One particular example we have 
considered is the preparation of an Einstein-Podolsky- 
Rosen (EPR) state in the positions and momenta of a 
pair of distantly-separated trapped atoms ||l3|] , which in 
turn leads to a scheme for the teleportation of motional 
quantum states fl^ . 

For the proposals presented in entanglement of 

the atoms' motional states is achieved through the trans- 
fer of entanglement from the quantum-correlated output 
light fields from a nondegenerate parametric amplifier. 
However, as pointed out in the trapped-atom cavity 
QED setup can itself act as a source of nonclassical light 



fields; nonclassical motional states (prepared by some 
means independent of the atom-cavity coupling) can be 
transferred to the propagating output field from the cav- 
ity, which can in turn be coupled to the motion of a 
second atom confined inside another, distant cavity. The 
potential of such a trapped-atom nonclassical light source 
is underscored by the impressive control with which the 
motional quantum states of single trapped atoms can be 
controlled and, more especially, by the variety of non- 
classical motional states that have in fact been prepared 
experimentally p^-p^. 

The purpose of the present paper, then, is to expand 
upon this idea and present schemes whereby "local" en- 
tanglement of orthogonal motional modes of a single 
atom (or of a few atoms) trapped inside an optical cavity 
is transformed, via propagating light fields, into "non- 
local" entanglement of the motional modes of distantly- 
separated atoms. In this way, the need for nonlinear 
(quantum) optical devices, such as nondegenerate para- 
metric amplifiers, is eliminated, and all of the desired 
operations are achieved using only trapped-atom cavity 
QED configuratons. 

In Section II we present a relatively simple example 
of a scheme for the manipulation of the motional state 
of a trapped atom (or atoms) in two dimensions using 
only (external) laser fields. This scheme allows entangle- 
ment to be established between two orthogonal modes 
and the results of some numerical simulations for the 
case of a two-mode "squeezed" state of the motion are 
presented. Section III describes the trapped-atom cavity 
QED system that facilitates the motion-light coupling 
and is, of course, central to our proposals; through this 
coupling, the quantum state (i.e., the entanglement) of 
one of the two relevant motional modes can be "mapped 
onto" a propagating light field and transferred to a dis- 
tant site. We examine this system in some detail, pre- 
senting numerical calculations that further the original 
analysis of looking more closely, for example, at the 
Lamb-Dicke (tight confinement) assumption involved in 
the model. Having described means for preparing "lo- 
cal" entanglement of motional modes and for transfer- 
ring motional quantum states between distant atoms, the 
preparation of distributed "nonlocal" motional state en- 
tanglement follows naturally, and in Section IV we dis- 
cuss several examples, such as the position-momentum 
EPR state and delocalized mesoscopic states of the mo- 
tion. 
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II. ENTANGLING ORTHOGONAL MOTIONAL 
MODES AT ONE LOCATION 

Schemes for coupling and manipulating orthogonal mo- 
tional modes in a single-atom trap have been put for- 
ward by a variety of authors p^-p6[| , typically involving 
stimulated-Raman transitions between internal atomic 
states. Here, we describe one such scheme, outlining how 
certain effective interaction Hamiltonians (of particular 
interest to us) may be realized and presenting some nu- 
merical simulations. This scheme involves only a single 
internal atomic level (after adiabatic elimination of the 
excited internal state) and generalizes to two dimensions 
a technique already used in ion-trap experiments to pre- 
pare, for example, one-dimensional squeezed states of the 
motion 1 15|. 



A. Coupling orthogonal motional modes of a single 
trapped atom 

In this scheme, counterpropagating laser beams 
aligned (in our case) in the x-z plane and coupled to 
the same internal atomic transition induce a coupling be- 
tween orthogonal {x and z) motional modes which, de- 
pending on the detuning between the two laser fields, 
may take the form of a linear mixer or a nondegenerate 
parametric amplifier. 

The physical setup is depicted in Fig. 1 (a) . Mathemat- 
ically, the situation is described by the Hamiltonian 

H{t) = ^AoiCT+(T_ -I- ^ fivj{b]hj + 1/2) 

j=x,z 

+ ih[El{x,z,t)a^ - Ei^{x,z,t)a+\ , (1) 

where Vx and are the harmonic oscillation frequencies 
along the x and z axes of the trap, hj are annihilation 
operators for the quantized atomic motion, and ct- = 
\g){e\ is the atomic lowering operator for the \g) ^ \e) 
transition. The total field incident on the atom is given 
by 

Ei^{x, z, t) — Ei{x, z, t) + E2{x, z, t) 



_£_ 

71 



(2) 



where a and P are determined by the angle of incidence 
in the x-z plane and satisfy (a^ + (3'^)^^^ — 1. The po- 
sition operators are given hy x = [fi/{2mi>x)Y^'^{bx + il) 
and z = [h/{2mv^)Y/'^{b^ + h\), and k = 27r/A is the 
wavenumber (taken to be the same for both fields). The 
detunings Aoi and 821 are given by Aoi =0^0—^1 and 
821 = LU2 — LUi, where luq is the atomic transition fre- 
quency and uji and UJ2 are the frequencies of the two 
counterpropagating fields. 

With the assumption that Aqi is large, the atomic ex- 
cited state can be adiabatically eliminated and sponta- 
neous emission neglected. In the Heisenberg equations of 
motion one makes the substitution 



.EL{x,z,t) 
a-^i , 

which gives, for the motional mode operators 
2i 



^01 



\Ei^{x,z,t)Y 



+ - 



-2ik{a.x-\-f3z)-\-i62it—i(p 



h.c. 



(3) 



(4) 



(5) 



where rj!^ — arjx and 77^ — f3riz, with rj^ — k{h/2mvx)'^^'^ 
and r^z — k{fi/2mVzY^'^ ■ 

Moving to a rotating frame with respect to the motion, 
i.e., defining bx.z = f'aj.xe""^"''^*, and making the Lamb- 
Dicke approximation, i.e., expanding the exponentials to 
first order in r/^ ^ , one derives 



-iS-2lt+i4- 



^01 

Aoi 
Aoi 



(and similarly for hz). This equation contains a variety 
of terms oscillating at a variety of different frequencies. 
With a judicious choice of the detuning 821 between the 
two laser fields, a particular term can have its time de- 
pendence removed, whereas the other (oscillating) terms 
can, provided the trap frequencies are sufficiently large 
(and different), be neglected in a rotating- wave approxi- 
mation. We now consider two cases of interest to us. 



1. Linear mixing: rotation of the motional state 

To realize a linear mixing of the x and z modes (i.e., 
a beamsplitter- type interaction), we choose the detuning 
between the two incident laser fields to be 821 = Vx~ Vz- 
With this choice, only one term in (^) is "stationary" , 
while all of the other terms retain oscillatory factors. 
Provided the frequencies Vx, Vz, and \i'x — Vz\ are large 
compared to the effective coupling strength between the 
two modes, these oscillatory terms can be dropped in a 
rotating-wave approximation, leaving 



Aoi 



2 „±lri 



(7) 



The effective interaction Hamiltonian is thus 
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with the interaction strength given by 



Aoi 



(8) 



(9) 



The physical process corresponding to this situation is 
illustrated in Fig. 1(b). Particular cases of interest are 
those in which the interaction is turned on for a time T 
such that = 7''/4 or xT = 7r/2 (we assume rectangular 
pulse shapes); the first case effects an operation equiva- 
lent to a 50/50 beamsplitter, while in the second case the 
result is a complete exchange of states between the two 
modes, i.e., for — n/2 one finds bx{T) — —hz{0) and 
h,{T)=K{0). 

Note also that, given suitable (unentangled) initial 
states of the x and z modes (prepared by some other 
means), one can generate entanglement with the beam- 
splitter {xT = 7''/4) operation. For example, an initial 
product of number states. 



1^(0)) = l^x). ® \N, 



(10) 



is transformed via a 50/50 beamsplitter operation into 
an entangled state of the form 



N 



(11) 



n=0 



where iV — + and {d„} are certain coefficients. 



2. Two-mode squeezing of the motional state 

Alternatively, we may choose 621 = v-x + i^z, h^ which 
case, following the same assumptions as above, one de- 
rives 



corresponding to 



'^V'xiz^'^e 



^01 



bl 



hx (fet&tei^-Hh.c), 



(12) 



(13) 



with X as before. The physical process corresponding to 
this situation is illustrated in Fig. 1(c); it is of course 
a parametric amplification process, which leads to "two- 
mode squeezing". Taking <j) — — 7r/2, the state produced 
after an interaction time T, given the modes are initially 
in their ground vibrational levels (achieved, for example, 
by sideband coofing), is 



IV'sq(r)) =5,,(r)|0),® |0), 
00 

= [cosh(r)]^"'^ [— tanh(7' 



where 



are Fock states of the motional modes and 



Sxzir) is the two-mode squeezing operator |27| , 



Sxzir) = exp 



(y 



bill 



(15) 



with r = xT- The entanglement between modes that 
is generated by this process is best expressed in terms of 
position and momentum variables. In particular, one can 
show for the state (^^ that the Wigner function in these 
variables is 

W{x,px] z,pz) 

= ^ exp{-[{i + Sf + {px~Pzf] e+2'-} 

X exp {- [(£ - zf + {px + p. f] e-^-} (16) 
^ C 5{x + z) S{px — Pz) for large r, (17) 



with C a constant |2^,g^. This, of course, corresponds 
to the original state considered by Einstein, Podolsky, 
and Rosen in their famous gedanken experiment pg| ] , al- 
though, in this particular situation the position and mo- 
mentum variables belong to the same particle. However, 
as we will describe below, through our scheme for motion- 
light coupling we are able to distribute this entanglement 
between two distantly-separated atoms. 



B. Numerical analysis 

Before describing the scheme for distribution of mo- 
tional state entanglement, we consider the results of some 
numerical simulations of the state manipulation schemes 
proposed above. By numerical simulation, we mean di- 
rect integration of the Schrodinger equation 



||V(i)) = ^i?ad(t)|^(0), 



(18) 



with iJad (t) given by (omitting the constant ground state 
vibrational energy) 



He,d{t) = huxblbx 



fivzb\bz 



(19) 



which is the form one obtains upon adiabatically elimi- 
nating the atomic excited state. Truncated number state 
bases are used to describe the harmonic oscillator modes 
describing the two-dimensional motion. 

We have focussed in our work on the two-mode squeez- 
ing, or EPR entanglement, operation of the previous 
subsection, determining in particular suitable parameter 
regimes for the efficient implementation of the scheme 
(requirements for the linear mixing scheme are basically 
the same). As a measure of how well the scheme works, 
we consider the fidelity 



F=|(V;.q(r)|V;(T)), 



(20) 



\m)z , (14) where |V(T)) 



is the simulated wave function after a 



m=0 



time r, and r = xT with x given by (^). 
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a. Results Sample values of F are given in Table I for 
a number of different combinations of trapping frequen- 
cies and Lamb-Dicke parameters. As one can see, two- 
mode squeezed states exhibiting significant degrees of en- 
tanglement (r > 1) can be prepared with high fidelity 
given sufficiently small Lamb-Dicke parameters and suf- 
ficiently large Vx, Vz, and \vx — Vz\ (compared to x)- 

It becomes computationally prohibitive for us to con- 
sider values of r much larger than 1.5, owing to the 
large number state bases required to cover adequately the 
population distribution of the two-mode squeezed state. 
However, the results presented here suggest that, for ex- 
ample, a state with r ~ 2 could be prepared with a fi- 
delity exceeding 0.9 for physically reasonable parameters. 

6. Experimental prospects If we take the example of 
trapped ^Be+ ions, rj'^ = v'z — 0.0707 and Vx/^z = 1/3 
correspond to actual trap frequencies v^j'^'K — W MHz 
and Vz/'i.-K = 34 MHz (which are essentially the frequen- 
cies achieved in the experiments of Monroe and coworkers 
|p6[), and thus, for the example given in the table, to a 
value of the interaction strength x/^tt = 22 kHz and a 
timescale for the preparation of the entangled state on 
the order of — MS. This should of course be much 
smaller than the timescale for motional state decoher- 
ence due to, for example, spontaneous emission, which 
we have neglected in our analysis. The rate of sponta- 
neous emission events infiuencing the motional dynamics 
(i.e., incoherent scattering) can be estimated to be on 
the order of (777^5^)/ Aqj^, where 7"^ is the atomic ex- 
cited state lifetime (see Appendix A, Part 1). Choosing 
7/A01 ^ 1, this rate can evidently be made much smaller 
than X- 



C. Two or more trapped ions: coupling collective 
and single ion modes 

Two or more ions confined in a linear ion trap in- 
teract strongly through Coulomb repulsion and their mo- 
tion along the axis of the trap is best described in terms 
of collective modes of vibration. Such collective modes 
have, of course, been fundamental to ion-trap quantum 
computer proposals |31|. In the present context, we wish 
to point out that the schemes of Section II. A could also be 
used to couple and entangle orthogonal collective modes 
of a linear chain of trapped ions. 

To illustrate this, we consider for simplicity just a pair 
of trapped ions confined in a harmonic potential and 
aligned along the z-axis, as illustrated in Fig. 2. As- 
suming very strong confinement in the x direction (so 
that the motion along this axis can be regarded as in- 
dependent of the z-axis collective motion) , the quantized 
(two-dimensional) motion of the ions can be described by 
the Hamiltonian 



Ha 



clz'coz 



hl^Rz 



'Rz 



CRz 



(21) 



where bj^, cqz, and crz are harmonic oscillator annihila- 
tion operators for single-ion motion along the x-axis, for 
the center-of-mass motion along the z-axis, and for the 
relative motion along the z-axis, respectively. The fre- 
quencies of the collective modes are related to the single- 
ion mode frequency by f oz = i^z and i/rz = V^i^z ■ 
To couple the various modes via laser light, excitation 
of only a single ion is necessary, and this excitation takes 
the same form as described by (|^) and (|2|); that is, we 
take 



^ion-laser— "^OlCr+ cr_ 



ih El{xi,zi,t)a'-_^^ -h.c. , (22) 



with 



Ei,{xi,zi,t) 



(23) 



Expressed in terms of the collective mode operators, one 
has M 



Akz 



(24) 



with the collective mode Lamb-Dicke parameters given 
by 



mz 



Vz 

V2' 



VRz 



Vz 



V2V3 



(25) 



As before, the atomic excited state is adiabatically elimi- 
nated and the position-dependent functions are expanded 
to first order in the various Lamb-Dicke parameters. The 
resulting Heisenberg equations of motion for the mode 
operators again display couplings between the modes 
with varying time dependencies, and again with a ju- 
dicious choice of the detuning 621 one can "select" a par- 
ticular coupling (neglecting the others in a rotating-wave 
approximation), provided the frequencies v^, vqz, and 
i^Rz are sufficiently large and different. The effective in- 
teraction Hamiltonians one can realize are thus (in the 
rotating frame) 



i?cff = < 





'blcoze^'f' - 


f h.c.) 


, S21 


= 


- VQz 


^Xo 


blxcle^^ - 


f h.c.j 


, S21 


= Vx 


+ V^z 


^XR 


(blcRze^"^ 


-1- h.c. 


) , S21 


= Vx 


- VRz 


^XR 




-t- h.c. 


) , ^21 


= Vx 


+ VRz 



with 



Xo 



^01 



XR = 



WxV'bJ^ 



^01 



(26) 



(27) 
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and 77^ = ar]^, 77^^ = pr]oz, v'rz = f^VRz- So, one again re- 
alizes linear mixers and parametric amplifiers, only now 
one of the modes is a collective mode for the two ions. 
With the addition of more ions, one introduces more col- 
lective modes, but the frequencies of these modes remain 
incommensurate and (for large i^z) well-separated 
implying that the above working could be readily gener- 
alized to three or more ions. 



III. COUPLING MOTION TO LIGHT 
A. Cavity-mediated light-motion coupling 

1. Model 

The basic setup we use to couple motion to light was 
originally considered by Zeng and Lin |Q. This setup 
consists of a two-level atom confined in a harmonic trap 
located inside an optical cavity. The atomic transition 
of frequency luq is coupled to a single mode of the cavity 
field of frequency uJcav and is also driven by an external 
(classical) laser field of frequency loa- The physical setup 
and excitation scheme are depicted in Fig. 3. The cavity 
is aligned along the cc-axis, while the laser field is incident 
from a direction along the y-axis (i.e., perpendicular to 
the X-axis). 

The Hamiltonian describing the atom-cavity system, 
including the atomic motion, takes the form (in a frame 
rotating at the laser frequency, uja) 

i?ac — ^ ^'^ji^lh + 1/2) + McAO^a -I- TjAoAfJ+CT- 
j=x,y,z 

+ n[EA{y,t)a+ + El{y,t)&_] 

+ %o sm{kx){a^a^ + (J+a) 

-fatte + tta + <7+ta + tt^_. (28) 



Here, {v.. 



X 1 1 ' 



y } are the harmonic oscillation frequencies 
along the principal axes of the trap, bj and a are annihila- 
tion operators for the quantized atomic motion and cav- 
ity field, respectively, = \g){e\ is the atomic lowering 
operator, and 5cA — ^cav — and Aqa = — (^a- The 
quantity EA{y,t) is the (possibly time-dependent) am- 
plitude of laser field A. The single-photon atom-cavity 
dipole coupling strength is given by g^, while the sine 
function describes the standing wave structure of the 
cavity field - we assume that the centre of the trap is 
located at a node of the cavity field. Finally, the last 
two terms in ( p8| ) describe the couplings of the cavity 
field mode and the atomic transition to "reservoirs" of 
external field modes (with Ta,c the "reservoir annihila- 
tion operators"), which produce damping of the cavity 
field and (free-space) atomic spontaneous emission, re- 
spectively Note that we will neglect any forms of 
motional decoherence associated with the trap itself. 



2. Elimination of the atomic excited state 

Heisenberg equations of motion are straightforwardly 
derived from the above Hamiltonian. Assuming the de- 
tunings of the light fields from the atomic transition 
frequency to be very large [i.e., Aqa ^ \EA\,go,5,Vj\, 
atomic spontaneous emission can be neglected and the in- 
ternal atomic dynamics can be adiabatically eliminated. 
In the equations of motion, this is done by making the 
replacement 



1 



iOA 



[EA{y,t) +go sin(fci)a] 



(29) 



in the equations describing the cavity and motional de- 
grees of freedom. The corresponding Hamiltonian then 
takes the form 

i^ac = X! T^Vj{h]bj + l/2) + h5^Aa^a 

i=x,y,z 

MEA{yM hgl . 2., ^^~t" 
^ sm [kx)a^a 



iOA 



iOA 



sm 

AoA 



(kx) [EAiy,t)a^ + El{y,t)a] 



(30) 



3. Lamb-Dicke and rotating-wave approximations 

The size of the harmonic trap is assumed to be small 
compared to the optical wavelength (Lamb-Dicke ap- 
proximation); this enables the approximations sin(fc2;) ~ 
77,(6, +6t), and EA{y,t) ~ EA{t)e-"l"' [with £a(0 real]. 
This second approximation would follow, for example, if 
the laser field forms a standing wave with the trap cen- 
tered at an antinode, i.e., with i?A(y) oc cos(fcy) ~ 1, for 
77y ^ 1 (for further discussion of this approximation, see 
Appendix B). 

To second order in 77,, quantum Langevin equations 
for the field and x-dimension motional modes are then 



-(k + iScA)a + i 



7 2 
Tx% 

AoA 



{bx + hl) 



bx = -ivxbx + i- . 

Aqa 

-iVL(t) (ate"^"^* 
where we have defined 

Vtit) = 



(bx + bi) 



ae 



'ainit), (31) 



(32) 



Vxgo£A{t) 



iOA 



(33) 



while the operator ain(0 obeys the commutation relation 
[ai„(t), aj'jj(i')] = S{t—t') and describes the quantum noise 
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input to the cavity field from the external field (in a frame 
rotating at the cavity frequency). The parameter k is the 
cavity field decay rate. 

Next, we choose the detuning between the cavity and 
laser fields to be ScA = Wcav — wa = i^x- In the above 
equations this results in "resonant" and "non-resonant" 
terms. With the assumption that the trap frequency Vx is 
large (which is consistent with the Lanib-Dicke assump- 

tion, since rj^ oc Vx ), such that i/^ 2> k, |r2(t)|, the 
non-resonant, or counter- rotating terms can be neglected 
in a rotating-wavc approximation. This leads to the pair 
of equations 

, , . .2r]lgl 1\ „ 

a = -(k + lUxja + I W^bx + 7: ] a 

^OA V 2/ 

^ in{t)e-'">"^k - V2^e-"'''*ai4t), (34) 
bx = -ii^xk + i^^^^jr^a^ahx - in{t)e"^''a. (35) 

The terms of second order in r]^ describe phonon- or 
photon-number-dependent phase shifts; these will in gen- 
eral be very small and can be neglectedf although they 
are retained for numerical calculations) |35| ] , which means 
that the effective interaction between the cavity and mo- 
tional modes is simply a linear coupling of the form 

H^s = hn{t) (atfo^e-*"^* + biae"*"'^ . (36) 



4- Adiabatic elimination of the cavity mode 

Although not essential for our purposes, a further sim- 
plification of the dynamics is possible if the decay rate 
K of the cavity field is sufficiently large that the cavity 
mode dynamics can also be adiabatically eliminated. In 
particular, if k ^ 1^(01 (but still with i/^ ^ k), then the 
equation for the motional mode reduces to 

hx ^ -[T{t)+tJ^x]bx + e'^^y^2me-'''^'ain{t) , (37) 
where we define 

T(t) = n{t)'^/K. (38) 

The motional dynamics thus reduces to that of a simple 
damped harmonic oscillator, with the (possibly) time- 
dependent damping rate r(t|. However, the quantum 
noise operator appearing in ( |37| ) corresponds to the light 
field incident upon the cavity, and hence one realizes a 
simple coupling of the motional mode to propagating light 
fields external to the cavity. More precisely, from the 
input-output theory of optical cavities ]27| , |36|| , it can be 
shown that the cavity output field is given, under the 
present circumstances, by 

aout{t)^-ain{t)- ^2T{t) 6,(i) . (39) 



where = e^^'^^bx^ and we have set (/)a = for sim- 
plicity. So, given a vacuum input field to the cavity, the 
output fight field is determined by the motional state of 
the atom confined inside the cavity. In this way, nonclas- 
sical motional states can be converted into nonclassical 
light fields; for example, entanglement between the a;- 
dimcnsion motional mode and, say, the z-dimension mo- 
tional mode can be converted into entanglement between 
the z-dimension motional mode and the propagating out- 
put light field. This light field may then be coupled to 
another atom-cavity system to generate distributed mo- 
tional state entanglement. 

B. Motional state transfer between distant locations 

Following work by Cirac et al. ||37[| on the transmis- 
sion of a qubit between two nodes of a quantum net- 
work, it is shown in that if the output field from one 
of our atom-cavity configurations is incident on a sec- 
ond such atom-cavity configuration, with the coupling 
between systems being unidirectional^ then with suitably 
tailored laser pulses £ai (t) and £a2 [t) applied at the two 
sites (amounting essentially to impedance matching, such 
that all of the light exiting the first cavity is absorbed by 
the second cavity) one may realize the motional state 
transfer 

Wx^ ® Wx^ - |0)i^) ® |</>)i'^ , (40) 

where \4>)x is an arbitrary quantum state describing the 
motion along the x-axis. Note that the cavity fields begin 
and end the transfer in the vacuum state. 

The scheme is able to operate in the regime where Q,{t) 
and K are comparable (in which case the transfer rate 
is essentially determined by k), but again the analysis 
simplifies in the case where the cavity modes can be adi- 
abatically eliminated as above. In this case, a master 
equation for the reduced density matrix describing the 
motional states of the two atoms, p,j., can be derived (in 
the rotating frame) as ||l^ 

Px = ri(*) {2bixPxb\x - b\xbixPx - Pxb\xbix^ 

+T2{t) {2b2xPxb\x - b\J>2xPx - Pxb\xb2x^ 

+ [PxbWb2xV^*^''*^-^] . (41) 

Example pulse shapes for this regime, specified 
through the effective coupling rates of the motional 
modes to the external light fields, Viit) and r2(t), are 
(taking 0ai = 4>A2) Q 

riW-r ^^^^"_^^ , T2{t) = T,{~t)^ (42) 
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assuming the transfer starts at t = — oo and concludes 
at t = +00, with r a constant. Armed with this ca- 
pabiHty, we are able to distribute quantum states of a 
material oscillator, and generate entanglement, between 
macroscopically-separated locations. 



C. Numerical analysis 

1. Damped harmonic oscillator model 

The description of the motional mode dynamics in 
terms of a linearly-damped harmonic oscillator coupled 
to propagating light fields. Eg. , represents a tremen- 
dous simplification of the model and offers a very direct 
and transparent scheme for state transfer between mo- 
tion and light. It is important then to gauge the validity 
of this simplification; here we present some results from a 
numerical analysis of the model starting from the Hamil- 
tonian 



i?ac= ff^x [blbx + 1/2 j + hS^Aa^a 



AoA 



A 



OA 



AoA 

and master equation 



sin(fci) (e-^^*at-he^^*a). 



(43) 



-[iJac! /o] + {2apa) — a) ap — pa) a) , 



(44) 



where the last term in 
The Hamiltonian (||) is t 



describes cavity dampin 
le one-dimensional form of (pi 
and neglects the position dependence of (based on ear 
lier arguments; see Appendix B). These equations make 
no assumption about the Lamb-Dicke parameter in the 
X direction and retain the dynamics of the cavity mode. 

We solve the master equation ( p^ numerically using 
truncated number state bases for the (harmonic oscilla- 
tor) cavity and motional modes. The particular example 
we have concentrated on is the decay of an initial coher- 
ent state of the motional mode, since this allows a simple 
analytical solution in the ideal case where the motional 
mode dynamics is exactly described by (^). In particu- 
lar, for an initial coherent amplitude a, the state evolves 
as \il){t))a: = e-*''-*/2|ae-(^''-+^)*)^ in the ideal case. 

A natural quantity to consider then in comparing the 
simulated evolution with the ideal behavior is the fidelity 



f{t) = 



(ae-(-^+r)*|p^|ae-(^''^+r)*) 



(45) 



where Px = Trcav{/5} is the reduced density operator for 
the motional mode and 



K \ AoA 



(46) 



a. Results Results for two values of the Lamb-Dicke 
parameter, rj^ — 0.1 and r]^ — 0.15, are shown in Fig. 4, 
where we plot the magnitudes of the motional mode and 
cavity mode amplitudes as a function of time, with the 
coherent state amplitude chosen to be a = VTO. In di- 
mensionless units we choose k = 1, I'x — ScA — 10, 
.gg/AoA = 0.2, and ry^jgo^A/AoA = 0.1 (< k), corre- 
sponding to r = 0.01. In Fig. 4(a), \{bxit))\ is seen to 
follow very closely the ideal behavior of a decaying coher- 
ent state, while Figs. 4(b,c) demonstrate that the cavity 
mode, after an initial transient period, adiabatically fol- 
lows the motional mode. There is a small but noticeable 
improvement (with respect to the ideal behavior) with 
the decrease in rjx from 0.15 to 0.1. Further improve- 
ments occur with smaller (rj^goSA/^aA)/!^ (condition of 
adiabaticity) and larger v^/k (rotating-wave approxima- 
tion with respect to the trap frequency), as one would 
expect. 

The fidelity f{t) is plotted in Fig. 5, from which it is 
clear that, for 77 = 0.1, the motional state remains close 
to the desired state at all times (at large t, of course, 
the states all approach the ground state). For rjj; = 0.15 
the deviation is more significant and the mapping of the 
motional state onto the light field is evidently degraded. 

b. Lamb-Dicke approximation Let us return briefly to 
some of the approximations made in deriving the ideal 
model. If we consider the Hamiltonian (^), the term 
proportional to Jq/Aoa is essentially negligible under 
the present circumstances due to the smallness of rj^ and 
of the intracavity photon number ((a^a) < 0.1 at early 
times and decreases with time). Consider now the ex- 
pansion of sin(fca:). To third order in rj^ 



I) sin(fci) = (S, + bl)-^ {b=c + bl) 

2 



Vx 



2 

2 2^^- " 6 ^ 



(47) 



The reduction of the atom-cavity dynamics to a coupling 
of the form a^bx + h.c. requires that the terms in the 
brackets proportional to ry^ have negligible effect (com- 
pared to 1). The combination of small rjx and rapid ro- 
tation means that the terms Ti^b^/6 and rj^Wx 1^ should 
be negligible, leaving a condition of the form 



\vl{^ + rix) « 1, 



(48) 



where fix = {b\.bx). This condition is reasonably well 
satisfied for the numerical examples considered, but a 
stronger condition would take into account the width of 
the number state distribution, i.e., the fact that vibra- 
tional population can reside in appreciable amounts in 
number states \n)x such that n > fix- If crfi is the 
variance of the number state distribution, then such a 
condition might take the form 



Tlx + aanJ < 1, 



(49) 
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where a 2 — S (i.e., several standard deviations from the 
mean). For the coherent state considered above, (7n_^ = 
fn^ = a/To, leading to a condition (taking a = 3) 

1 



3a,-,J ~ 1077^ « 1. 



(50) 



With = 0.15, lOr/^ = 0.225, and so one might ex- 
pect that the Lamb-Dicke approximation starts to break 
down. 

In the case of a two-mode squeezed state of the motion, 
the number state distribution of the individual modes is 
akin to that of a thermal mode, for which (Jn^ = [n^ + 
n^Y^"^ — "a: + 1/2 for n.j. > 1. With a = 3, we then 
require that 

1 



-ry^(4n, + 5/2)«l. 



(51) 



which, for r\x — 0.1, reduces to — sinh^(r) ^ 50. This 
condition is reasonably well satisfied for values of r up 
to ~ 1.6; beyond this, a smaller value of r]^ would be 
desirable. 



2. Motional state transfer 

The next problem we wish to examine in more detail is 
that of transferring a motional quantum state from one 
atom to another at a distant site using the procedure 
outlined in Section III.B and depicted schematically in 
Fig. 6. The coupling laser fields, and hence the equa- 
tions of motion, now have an explicit time dependence. 
This, combined with the increased dimensionality of the 
problem (now that wc have two atom-cavity systems to 
simulate), leads us to employ the technique of Monte- 
Carlo wave function simulation (see, for example, [^). 

In fact, the Monte-Carlo wave function technique is not 
only convenient, but also highly appropriate for model- 
ing and analyzing the state transfer procedure as shown 
in Fig. 6. In this approach, the evolution of our quan- 
tum system is simulated by propagating a wave function 
|\['(i)) according to the Schrodinger equation 



d_ 

di 



-H,sit)\^it)), 



(52) 



where Hesit) is a non-Hermitian effective Hamiltonian. 
For the state-transfer configuration we are considering 
here, our effective Hamiltonian takes the form 

— i'hKia\ai — i'hK.2d\d2 — 2ifi^ kiK2 d\ai, (53) 

with 



^9oi . 



A 



OA 



sin(fc%) (e-^^-^a] +e^*-^a,) . (54) 



[For simplicity, we now omit constant energy shifts from 

the Hamiltonians Hii' (t) , and the frequencies of the two 
coupling lasers are assumed to be the same (wa).] The 
second line in (|5^) follows from the cascaded-systems for- 
malism in which one assumes a unidirectional coupling 
between the cavities |^,|8|,|9| . 

The propagation described by ( |52| ) is interrupted at 
random times {tr} by wave function collapses, or quan- 
tum jumps. 



(55) 



(which include renormalization of the wave function), 
where C is an appropriate jump operator and the proba- 
bility density for a jump to occur during the time interval 
from t to t + dt is {^>{t)\CW\'^{t))dt. For our situation 
the jump operator is 



(56) 



which can be identified with the destructive detection 
of a photon by the photodetector monitoring the output 
channel from the second cavity. 

For ideal quantum transmission we require that a 
quantum jump (that is, a photon detection) never oc- 
curs, i.e., C\'9{t)) = for all times t (which also means 
that the norm of the wave function remains equal to 1 
at all times). In other words, all of the quantum infor- 
mation is transferred from atom 1 to atom 2 and none 
is lost to light fields propagating away from the system. 
The laser pulse profiles derived from the simplified mas- 
ter equation model (|4l| ) should approximately satisfy this 
condition and facilitate high-fidelity state transfer, pro- 
vided the various parameters of the system are chosen 
appropriately. Here we want to assess more carefully the 
performance of these profiles for a more comprehensive 
model of the atom-cavity dynamics, as described by the 
Hamiltonians in (^J). 

We assume an initial state of the form 



|*(t = ~oo)) 



|0)W 0|O)E®|O)i^\ (57) 



where \(j))x is the particular motional quantum state to 
be transmitted. The target state of the transmission is 
then 



\^{t = +00))targct - |0)il) ® |0)« (g> |0) 



(58) 



For simplicity, we will also assume identical atom-cavity 
systems, i.e.. 



1^x1 = 1^x2 = I'x, Vxl = r]x2 = TJx, 
gOl = 502 = 50, Kl = K2 = K, S, 



(1) _ x(2) 



3c A. 



(59) 



8 



a. Laser pulse profiles The time dependence of the 
laser field £ai {t) is chosen to satisfy 



A, 



with 



OA 



r = 



ri(t) = r- 



max 
Al 



iOA 



(60) 



(61) 



while £A2{t) is chosen such that r2(i) = ri{—t). Note 
again that these example forms for the temporal profiles 
of the laser fields are derived from a theoretical analysis 
in which the dynamics of the cavity modes are adiabati- 
cally eliminated, and hence, depending on the choices of 
parameters, they may not be optimal choices. We will 
return to this point when we present the numerical re- 
sults. 

b. Truncation of trigonometric functions To help 
speed up our computations we have in fact performed 
most of our simulations using expansions of the opera- 
tors sin^ (kxj) and sm{kxj) that are truncated to third 
order in the Lamb-Dicke parameter rj^. For the range of 
parameters considered here, we find little difference be- 
tween results obtained using the full trigonometric forms 
and those using the truncated versions. With the trun- 
cated versions it is also straightforward to move to a ro- 
tating frame so as to remove the systematic (fast) evolu- 
tion given by the first line in (|^. This does of course 
result in rapidly rotating factors multiplying some of the 
remaining terms of the Hamiltonians [in particular, fac- 
tors of the forms e±'(''-+'^=A)* or e±*(3i^.±5cA)t]^ ^nd we 
retain these terms in the simulations. 

c. Fidelity of the transmission In the context of 
our wave function simulations, the ideal situation cor- 
responds to the case where, following the laser pulse se- 
quence, the norm of the wave function is still equal to 
1, and hence no jumps have occurred. However, due to 
non-ideal operating conditions (for example, rapidly ro- 
tating "off-resonant" terms making a finite contribution 
to the dynamics) the norm of the wave function does de- 
cay and hence there is a finite possibility of a photon 
detection and consequent loss of information. Now, in 
a typical application of the Monte-Carlo wave function 
approach one averages over many trajectories to obtain 
a density operator for the system. From this density op- 
erator one could, in the present context, compute the 
average fidelity for the quantum state transmission. 

However, the decay of the norm in a single trajectory 
already provides us with a good indicator of the perfor- 
mance of the transmission. If, for example, during the 
pulse sequence and state transfer the norm decays to a 
value of 0.9 (without any jumps occurring), then we can 
make the general statement that in 90% of our attempts 
the state will be transferred successfully. Note that, in 
the "no-jump" case, on renormalizing the simulated wave 
function |'I'(tf))sim (where <t is the finishing time for the 



simulation) we find that the fidelity of the transmission, 
defined to be 



F 



(vI'(tf)l*( + CX3))targetr, 



(62) 



is very close to 1 (i.e., > 0.99) for all of the numerical 
examples we consider below. The value of 0.9 then also 
essentially sets a lower bound on the average fidelity of 
the state transmission. In the 10% of cases where one or 
more quantum jumps do occur the final state |5'(tf ))sim 
may, depending on the precise nature of the state being 
transferred, still have a finite (and even substantial) over- 
lap with the target state, and so the average fidelity will 
actually be larger than 0.9. 

From a practical point of view, it is also worth not- 
ing that it would in principle be possible to post-select 
high-fidelity state transmissions by monitoring the out- 
put from the second cavity. If a photon is not detected in 
this output then the transmission is known to have been 
successful (assuming perfect detection efficiency). If a 
photon is detected then the system can be reset and the 
transfer attempted again (or some form of error correc- 
tion procedure could be applied). 

d. Results As the state to be transferred, we have 
concentrated on the following examples: (i) a truncated 
phase state. 



1 



N 



(63) 



i.e., a uniformly- weighted coherent superposition of the 
first + 1 Fock states, (ii) a pure Fock state 



(1) 



and (iii) a Schrodinger Cat state 



1 



(64) 



(65) 



with \a)i^^ a coherent state and7V+ = [2(1 + e"2|"l')]i/2. 

In dimensionless units, we again set k = I and choose 
gg/AoA = 0.2 and igo£r"')/^oA = 1. For r;, = 0.1 this 
gives {T]xgo£A'^^)/^oA = 0.1k and so the adiabatic ap- 
proximation [{rj^goSA'^^)/ '^OA ^ k] used in deriving the 
laser pulse shapes should be reasonably good. Indeed, we 
find that the improvement in the fidelity of the transfer 
one obtains when rjx is reduced from 0.1 to 0.0707 (see ta- 
bles below) results primarily from an improvement in the 
Lamb-Dicke approximation rather than in the adiabatic 
approximation. 

(i) Truncated phase state: Some results illustrating the 
performance of the transfer for this state are presented in 
Tables H (TV = 10) and |l| (A^ = 20). We have considered 
various combinations of the trapping frequency fx and 
the Lamb-Dicke parameter rj^- The third column gives 
the norm of the wave function at the conclusion of the 
transfer operation (with quantum jumps "turned off"). 
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which, as discussed earlier, effectively amounts to a lower 
bound on the average transfer fideHty, i^avc- 

As one can see, a large trapping frequency is very im- 
portant for obtaining a high average transfer fidelity. In 
particular, a value Vx/ k, > 10 is necessary if one wishes 
to obtain a success rate exceeding ~ 80% for the case 
N = 20. As the dimensionality of the motional state 
being transferred increases (i.e., higher Fock states are 
populated), the size of rjx also becomes more critical. 
This can be seen by comparing the results in Tables || 
and [II. Smaller values of rj^ are clearly required (for a 



^20 



given value of Vx) in order for the state 1/ v 21 ^ 
to be transferred with a fidelity comparable to that for 
the state l/v^X^n'Lo l")^- '^^'^ relative improvement 
in performance with a decrease is r]x is also more pro- 
nounced for the state of higher dimensionality. 

With the inclusion of quantum jumps some interesting 
and complicated behavior is observed in the individual 
trajectories. For the phase state considered the overlap 
of the transmitted state with the ideal state can be sub- 
stantial (~ 0.40 — 0.75), but the precise nature of the 
transmitted state depends crucially on the time at which 
each jump occurs. 

(ii) Fock state: Results for transmission of the Fock 
state \n = 10)x are given in Table IV. The perfor- 



mance is in fact comparable to that for the phase state 
l/V2TX]n=o l"')^' which actually has the same mean 
phonon number of 10. 

The effect of quantum jumps on the transfer is, as one 
would expect, more severe than in the case of the phase 
state; overlap with the state \n = 10)x disappears com- 
pletely with a single photon detection. However, after 
the jump the resulting final state is not necessarily what 
one might naively expect (i.e., the state \n = 9)x)- It is 
possible, as a result of the finite effects of "non-resonant" 
terms in the dynamics, for the transferred state to take 
the form of a superposition of the states \n — 9)x and 

\n = n)x. 

(iii) Schrodinger Cat state: Finally, we consider the 
transfer of a mesoscopic superposition of coherent states 
of the form ^^:^^'^{\a)x + \ - a)x), with a = VTO. This 
state also has a mean excitation (phonon) number equal 
to 10 and the results shown in Table ^ are very similar 
to those of the previous two examples. 

The effect of quantum jumps is particularly interesting 
in this case, as we find that, with a single photon detec- 
tion, the transferred state is in general very close to the 
state A/r^(|a)a; — I — a)a;) (i-e., an "odd" Schrodinger Cat 
state), suggesting that an operation on the motional state 
amounting to an application of the annihilation operator 
bx2 would largely restore the original state. Note that the 
even and odd Schrodinger Cat states J^^^{\a)x + \ —cy)x) 
and MZ^{\a)x — \ — oi)x) have been proposed, in the con- 
text of quantum computation, as logical qubit encodings 
for the correction of bit-flip errors caused by amplitude 
damping (as occurs in our system) pO[. They have also 



been proposed for use in a secure quantum key distribu- 
tion protocol Q . 

e. Atomic spontaneous emission While our model 
and simulations have included cavity damping, we have 
until now ignored any effects associated with atomic 
spontaneous emission, which arises from the small but 
finite probability for the atom to be in its excited inter- 
nal state |e). In Appendix A (Part 2) an approximate 
expression is derived for the rate at which atomic spon- 
taneous emission is expected to influence the motional 
state of the atom. In order for this rate to be much 
slower than the characteristic state transfer rate F of the 
present configuration, one requires that the condition 



K7 



> 1 



(66) 



be satisfied, where 7 is the atomic spontaneous emission 
rate. Hence, one desires the regime of strong-coupling 
cavity QED, for which g^/ln-y) > 1. 

/. Experimental prospects Let us focus again on the 
case of trapped ^Be~'' ions, although other ionic or atomic 
species could evidently be considered (and might perhaps 
be preferable from the point of view of cavity QED ex- 
periments). The ^Si/2 ^ ^Pi/2 transition wavelength is 
A = 313 nm and the transition linewidth is ^/{2tt) = 
19.4 MHz. If we assume, for example, that the mir- 
rors forming the cavity have radii of curvature equal to 
5 cm and are separated by a distance / = 1 (2) mm, then 
5o/(27r) = 5.3(3.1) MHz. For a cavity finesse of 75,000 
one obtains k,/{2tt) = 1.0 (0.5) MHz, and so 10.gg/(K7) = 
14 (10), while a trap frequency of VxHItt) — 22 MHz (cor- 
responding to a Lamb-Dicke parameter r]x — 0.1) gives 
Vxl K, = 22 (44). With these parameters one would antic- 
ipate a state transfer rate F/(27r) ~ 10 — 20 kHz. Note 
that the timescales for motional decoherence and heat- 
ing observed in recent trapped ion experiments are of the 
order of milliseconds or longer |l7|,|18l 



There are of course many other possible combinations 
of parameters which should satisfy (at least approxi- 
mately) the essential requirements of the scheme. Much 
larger cavity finesses may be possible, but challenging, 
at this wavelength, allowing smaller mirror separation 
and larger dipole coupling strength go (which increases 
as the cavity mode volume decreases). However, it may 
be difficult in practice to bring mirrors very close to- 
gether about an ion trap (for example, because of the 
technical problem of charge build-up on the mirrors). 
For this reason, it is perhaps advantageous to consider 
neutral atom experiments (in particular, alkali atom ex- 
periments), where, in fact, optical dipole traps (see, for 
example, |43[| ) or microscopic magnetic traps (see, for 
example, p4|) should allow confinement of atoms in the 
Lamb-Dicke regime with trapping frequencies also in the 
MHz range. Furthermore, cavity QED experiments with 
alkali atoms are already highly developed (for recent ex- 
periments, see 1 43 4^,^) and spectacular cavity finesses 
can be achieved at the relevant (longer) wavelengths 
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D. Multiple cavity-confined atoms 

It is interesting to consider the possibility of having 
more than one atom coupled to the cavity field at a 
time. The kind of situation one might imagine is like 
that realized in the recent experiment of Ye et al. , in 
which (single) atoms were trapped inside a microscopic 
optical cavity using a far-off-resonance dipole-force trap 
(FORT). The FORT was actually produced by excitation 
of an "auxiliary" longitudinal cavity mode, and so had 
a standing-wave structure colinear with (but of slightly 
different periodicity to) the cavity QED mode of interest. 
Hence, given sufficiently strong confinement, the FORT 
can be regarded as a chain of individual cavity-confined 
microtraps, each of which could in principle be occupied 
by a different atom. 

With, for simplicity, a single coupling field £a incident 
on the atoms from the side of the cavity, a Hamiltonian 
for the system can be written as (assuming a single micro- 
trap frequency i^x , and tight confinement along transverse 
directions as well) 

Ha.c= ^ Tt-Vx (bljbxj + 1/2^ + hScAO.'^a 
j 

_ y M sin2(fcx, + 0,)d^a - T ^ 



E 



■ sm 



iOA 



{kxj + e.j) {e-"l"'a) + e'-^^a) , (67) 



where 9j defines the position of the center of the j-ih 
microtrap relative to a node of the cavity QED field. 
In the Lamb-Dicke limit we can write 

siTi{kxj + 9j) = sin(fcxj) cos{9j) + cos{kxj) sin(6'j) 

^ Vx (bxj + blj^ cos(0j) + smidj). (68) 

Assuming that go/^OA <t:. 6 — and that go{a^a)^^^ <C 
£a, we can neglect the term in ( |67| ) proportional to 
^q/Aqa, while in the last term we can neglect the con- 
tribution proportional to sm{9j){e~'^'^^a^ + e^'^-^a) as this 
term is rapidly rotating (at frequencies ±^cA, with 6cA = 
Ux ^ 9o£a/ ^oa) compared to terms of the form a^bxj 
and hi., a. 

Consequently, we can describe the dynamics of the sys- 
tem through equations of the form 



a = —{k + iScA)a 



-le 



2Ke 



bxj = -ivxbxj - ie^^^^jCL, 



'aUt). (69) 



(70) 



with 



(71) 



Adiabatic elimination of the cavity mode as before then 
leads to 

bxj = -(Fj -I- ivx)bxj - ^ ^^^^ bxk 

+e^'^*y2l7e-^''-*ain(0, (72) 



where F^ = H'j/k. So, through the laser and cavity fields 
one realizes a coupling between the motional modes of 
different atoms. 

Now, defining the "collective" mode operator 




cos{9k)bxk, 



(73) 



with Ncs = X]j COS^(0j) (such that [Bx,Bl] = 1), one 
finds 



Bx = -{N,sT + iVx)Bx - V2iWe-'"-*ain(i), (74) 

with F = {rjxgo^A/ ^oa)'^ / K. This points to the possibil- 
ity of using the collective mode for quantum state trans- 
fer and storage Q. A potential advantage of such an 
approach is the enhanced effective transfer rate A^offF, 
which (for iVoff > 1) should enable the requirement of 
strong-coupling cavity QED to be relaxed somewhat. 

The possibility also arises for the preparation of some 
very interesting many-atom entangled states. For exam- 
ple, motional state transfer of a Schrodinger Cat state 
N'+^i\a)x + \—a)x) to the collective mode, perhaps from a 
single-atom (cavity QED) source, would result in a many- 
atom state of the form (omitting the normalization) 

\ai,a2,a3, ■ ■ ■)x + \- ^i, -"2, -"3, ■ • ■)x, (75) 



where aj = acos{9j)/y/NeS. 



E. Coupling to collective trapped ion modes 

Earlier, we discussed a scheme for coupling orthogo- 
nal single-ion and many-ion collective vibrational modes 
using a pair of laser fields. This scheme could also be uti- 
lized in the context of state transfer and storage, albeit 
somewhat indirectly, as depicted in Fig. 7. In partic- 
ular, state transfer between motion and light and from 
one cavity location to another could be achieved using a 
single-ion interaction with the field mode (assuming that 
individual addressing of the ion by the coupling laser is 
possible) , after which the ion-cavity coupling is switched 
off and the coupling between the single-ion mode and a 
particular collective motional mode is switched on (for 
example, to transfer a quantum state from the single-ion 
mode to the collective mode) . Storage of quantum infor- 
mation in certain collective vibrational modes could be 
advantageous due to the very slow heating rates experi- 
enced by these modes [Eoi. 
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IV. ENTANGLING DISTANT ATOMS 

Having described and analysed basic techniques for 
preparing "local" entanglement of vibrational modes, and 
then for coupling vibrational modes to propagating light 
fields, we now want to outline some specific and interest- 
ing possibilities for distribution of entanglement between 
distant locations. The procedure is straightforward: one 
prepares an entangled state of the orthogonal vibrational 
modes in the x- and, say, z-directions, after which the x 
mode is coupled to the cavity and thence to the external 
light field. With appropriate time-dependent coupling 
we have seen that the properties of the x mode can be 
transferred with high fidelity to the corresponding mode 
of a distant cavity-confined atom (or atoms). That is, we 
can perform the transformation 



where {z^^\p^P} and {x^'^\p^'} are position and mo- 
mentum variables for the atoms at locations 1 and 2, 
respectively. 

So, one prepares an EPR state in position and mo- 
mentum of a pair of (distantly) separated atoms. Apart 
from being of great historical significance, this state also 
constitutes the essential resource for the teleportation 
of continuous variables 51 1. In particular, in the 

present context it offers the possibility of teleporting 
atomic center-of-mass wave functions between distant 
sites llTl. 




(76) 



where \n)x\ is a vibrational Fock state of the atom at 
location i. Hence, the x and z modes of the atoms at 
locations 2 and 1, respectively, become entangled. 

Note that theoretical proposals exist for schemes that 
would in principle allow the preparation of arbitrary 
two-dimensional vibrational states of a trapped atom 
p2|- p4| , p6[ . It follows that, using the state transfer 
scheme, it would in principle be possible to prepare an ar- 
bitrary entangled state of distantly-separated atoms. Be- 
low we highlight briefly just a few examples of particular 
interest, with emphasis on possibilities associated with 
the particular two-mode entangling schemes presented in 
Section II. 



A. Examples 

1. EPR state preparation between distant atoms 

Here, we use the scheme of Section II. A. 2 to prepare 
the state |-0sq(?'))*-^^ given in (p^, following which the 
state transfer operation outlined above produces the de- 
localized state 



cosh ^(r) [— tanh( 



(77) 



The corresponding Wigner function is 
= -i exp {- [(i(2) + z(i))2 + [j^ - p^yf] 

X exp {- [(5(2) _ |(1))2 + (p(2) + p(l))2] g-2.r| ^ (73) 



2. Delocalized mesoscopic states 

As mentioned earlier, the linear mixing operation de- 
scribed by (^ can also generate entanglement given suit- 
able initial states. Consider, for example, the situation in 
which, say, the z mode is initially prepared in its ground 
state while the x mode is prepared in a Schrodinger Cat 
state, i.e., 

|V'(0)) = -^(|a)«-H|-a)W)®|0)«. (79) 

Application of the linear mixing operation for a time T = 
7r/(4x) will transform this into the state 



1 



a 

72 



(1) 



—a 
V2 



(1) 



—a 



(1) 



a 



(80) 



The state transfer procedure then generates the delocal- 
ized state 





\ (2) 




\ (2) 


a 


( 


) ® 


V2A ^ 




V2 



(81) 

From what we have seen in our numerical analysis, the 
coherent state amplitude a can be reasonably large (while 
still allowing high fidelity operations and transfers), and 
so the entangled state ( plj) could be regarded as a delo- 
calized mesoscopic state [|52| . 

Schemes have also been proposed for preparing a state 
of the form 



1 
72 



where \N)'^x'^ is a Fock state (see, for example, 
would lead to a delocalized state 



1 
71 



|iv)(2) « |o)(i) -f |o)(2) ® |7V)ii: 



(82) 
, which 

(83) 
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Such states are of potential interest in the context of 
phase sensitivity in a two-mode interferometer, where 
they should allow measurements at the Heisenberg un- 
certainty limit (see, for example, and references 
therein) . 



entangled and delocalized states of the form considered 
here would offer some unique and fascinating opportuni- 
ties for fundamental tests of quantum mechanics versus 
local realism using massive particles |13|,|52|,|61|. 



3. Many-atom entangled states 

Through the many-atom configurations discussed in 
Sections III.D and III.E one can generalize the entan- 
glement distribution procedure to collections of trapped 
atoms or ions located at the distant sites. For example, 
one could imagine entangling collective vibrational modes 
of two strings of trapped ions at separate locations. Such 
schemes complement proposals for entangling collective 
internal atomic states of separated free-space atomic en- 
sembles using propagating light fields (see, for example, 
1). 
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APPENDIX A: EFFECTS OF ATOMIC 
SPONTANEOUS EMISSION 



V. CONCLUSION 



1. Motional state preparation 



In this work we have described and analysed schemes 
(i) for producing motional state entanglement of an atom 
(or atoms) at one location and (ii) for distributing this 
entanglement between atoms at distant locations. The 
particular example we have focussed on for generating 
(local) entanglement was chosen for its relative simplic- 
ity and for its close relation to schemes already imple- 
mented in the laboratory. A further motivating factor 
which should also be emphasized is that the operations 
that it allows are of direct relevance to quantum com- 
munication and computing with continuous variables, for 
which squeezers and linear mixers are basic elements (see, 
for example, [|ll],|l]j57|-ra ) . 



The cavity-QED-based motional state transfer scheme 
that enables the entanglement to be distributed between 
distant sites has been introduced elsewhere ||l^, but the 
present work extends significantly the numerical analy- 
sis of this scheme. In particular, we have considered the 
transfer of what may be regarded as mesoscopic motional 
states and from this we have been able to gauge the valid- 
ity of some of the more fundamental assumptions implicit 
in the scheme, such as the Lamb-Dicke and rotating- wave 
approximations. Our calculations suggest that states of a 
substantial "size" can be transferred with high fidelity for 
physically reasonable parameter values and we have dis- 
cussed some possible experimental scenarios, encouraged 
by spectacular recent advances in atom and ion trapping 
technology and in cavity QED. 

The possibilities offered by the schemes for distributed 
entanglement are many and varied, and we have dis- 
cussed just a few examples in Section IV, along with some 
possible applications such as quantum teleportation of 
atomic wavepackets. Of course, in addition to potential 
applications in quantum communication and computing. 



Clearly, a very important assumption is that the effects 
of atomic spontaneous emission can be neglected. In a 
master equation approach, atomic spontaneous emission 
with the effects of recoil taken into account is modelled 
by a term of the form (considering, for simplicity, motion 
only along the x axis) [|62l 



{/5}spon = 2 (2'^-/5'3-+ - <3-+<T_/5 - pa+a-) , 



(Al) 



where 



-1-1 



P = - I du W{u)e' pe 



du I^(u)e*"''-(^-+^x)pe-™''-(^-+^i) 



(A2) 



Here, 7 is the spontaneous emission rate and W{u) — 
(3/4) (1 -I- u^) describes the angular distribution of spon- 
taneous emission for an atomic dipole transition. 

Staying in one dimension and considering the prepara- 
tion of a squeezed state of the motion, the appropriate 
combination of laser fields is 



Ei^(x,t) = Ei(x,t) E2{x,t) 



71 



, (A3) 



with 821 — Adiabatically eliminating the atomic 

excited state, which amounts to setting 



(J- 



.Ei^{x,t) 



^01 



(A4) 



in (Al), one finds that the leading order (in rf^) contribu- 
tion to the motional dynamics arising from atomic recoil 
due to spontaneous emission takes the form 



13 



where 



4 



1 - e 



-?i52ii |2 



II 



(A5) 



(A6) 



So, the rate at which atomic spontaneous emission can 
be expected to influence the motional dynamics is on 
the order of 77y^£^/AQj^, whereas the rate at which the 
squeezed state is prepared is essentiafly x = 4?7^5^/Aoi, 
leading to the condition 



An 



< 1 



(A7) 



for spontaneous emission to have negligible effect. 



2. Motion-light coupling 

The analysis of the effects of spontaneous emission in 
the context of the atom-cavity state transfer scheme fol- 
lows the above working, only now the relevant substitu- 
tion is 



.90 



sin(A:i;)d. 



(A8) 



Assuming that £a ^ ilxgo\/ (o^a), to leading order in rjx 
one obtains a term of the form 



lOA^A 



(A9) 



Hence, in order to be able to neglect the effects of spon- 
taneous emission on the transfer process, one requires 
that 



A ^ „2 7 



> Vx 



10 Ag^ 



K7 



>1. (AlO) 



This, not surprisingly, corresponds to the regime of 
strong coupling in cavity QED. 



APPENDIX B: LAMB-DICKE 
APPROXIMATION FOR THE ATOM-CAVITY 
COUPLING LASER 

In deriving our model for motion-light coupling we 
make the (simplifying) assumption that position depen- 
dence of the external coupling laser can be neglected, i.e., 
Ex{y-,t) ~ £x{t)e'^'^^. We mentioned that this could be 
justified, for instance, in the case where the laser field 



forms a standing wave with the trap centered at an antin- 
ode, i.e., Ej^{y,t) (x cos{ky) ~ 1, assuming tight confine- 
ment in the j/-direction. 

However, a traveling wave laser field (which is most 
likely a simpler proposition from an experimental point 
of view) should also suffice, as the following argument 
shows. Assume that 



EAiy,t) = £Ait)e 



iky 



(Bl) 



Then, first of aU, |i?A(y,i)P = £A{t)^- Secondly, on adi- 
abatically eliminating the atomic excited state, we have 
an interaction term of the form (omitting constant coef- 
ficients) 



^nyCby+bl)sin[7^.Abx + bl)]ia + a^). 



(B2) 



To second order in the Lamb-Dicke parameters, this ex- 
pression takes the form 



r]^{b^ + bl){a + a^) 

+iVx%{by + bl){b.^ + bl){a + a^) 



(B3) 



Now, we have 6cA — i^x, and assuming that i^x and Vy 
are both large, and that \2h'x ~ Vy \ is also large, the only 
contribution in (B3) that is not rapidly rotating is the 
desired term rjx{a)bx + b\.a). That is, the contribution 
from the position dependence of EA{y, t), as well as being 
of order 1]^%, is also rapidly rotating (assuming large i^y), 
and hence can be neglected. 
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TABLE I. Two-mode squeezed state preparation: exam- 
ples of fidelities achievable for several different trapping con- 
figurations. Note that, in the dimensionless units used here, 
the values of x correspond to the choice £^/Aoi — 0.1. 
Also, since rji oc Xj^/vl, setting = means that 

(a//3) = yju^/vz. 









X 


r ^XT 


F 


0.1 


1 


3 


0.004 


1 


0.991 


0.1 


1 


3 


0.004 


1.5 


0.932 


0.1 


1 


4 


0.004 


1.5 


0.955 


0.0707 


1 


3 


0.002 


1 


0.996 


0.0707 


1 


3 


0.002 


1.5 


0.975 


0.0707 


1 


4 


0.002 


1.5 


0.986 


0.0577 


1 


3 


0.00133 


1 


0.998 


0.0577 


1 


3 


0.00133 


1.5 


0.987 


0.0577 


1 


4 


0.00133 


1.5 


0.994 
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TABLE II. Numerical results for transfer of the trun- 
cated phase state l/'s/lTX]l=o l"')^:- Parameters are k = 1, 
ffg/AoA = 0.2, {go£T'^)/Aol = 1.0. 



__nx i^x{= Sca) (^ftf)l^(ff)) (no jump) 

0.1 5 0.65 

0.1 10 0.90 

0.1 20 0.96 

0.0707 5 0.66 

0.0707 10 0.91 

0.0707 20 0.97 



FIG. 1. (a) Coordinate axes and laser configuration for 
two-dimensional motional state manipulation of a trapped 
atom. (The trapping potential is not shown.) (b) Excita- 
tion scheme for the linear mixer interaction, (c) Excitation 
scheme for preparation of a two-mode squeezed state of the 
motion. Only the first few vibrational levels are shown, and 
the excited atomic state is omitted for simplicity. 

FIG. 2. Laser configuration for coupling two-ion collective 
and single-ion vibrational modes. 



TABLE III. Numerical results for transfer of the trun- 
cated phase state 1/V2l'^^"_^^ |n,),r. Parameters are k = 1, 
9^/X,x = 0.2. ((/„^_'r")/A„A = 1.0. 



Vx 


l^x{= Sca) 


(*(tf)|«'(tf)) (no jump) 


0.1 


10 


0.79 


0.1 


20 


0.88 


0.0707 


10 


0.84 


0.0707 


20 


0.94 


TABLE IV. 


Numerical results 


for transfer of the Fock 


state |n = 10) a:. Parameters are 


K = 1, go/AoA = 0.2, 


{goSri/^oA 


= 1.0. 




Vx 


l'x{= Sca) 


(*(tf)|*(ff)) (no jump) 


0.1 


10 


0.82 


0.1 


20 


0.92 


0.0707 


10 


0.85 


0.0707 


20 


0.95 



TABLE V. Numerical results for transfer of the 
Schrodinger Cat state N'^^Hajx + \ — o:)x), with a = vTo. 
Parameters are /t = 1, Qo/Aqa = 0.2, (gof a ^'')/Aoa = 1.0. 

Vx i^x{=Sca) (^(tf)l^(ff)) (no jump) 

0.1 10 0.81 

0.1 20 0.91 

0.0707 10 0.85 

0.0707 20 0.95 



FIG. 3. Schematic of proposed (a) experimental setup and 
(b) excitation scheme for state transfer between the motion 
of a trapped atom or ion and a quantized cavity mode of the 
electromagnetic field. All input and output to the cavity mode 
is through just one mirror, i.e., the other mirror is assumed 
to be perfect. 

FIG. 4. (a) Decay of the motional mode amplitude 
\{bx{t))\ with time, given an initial coherent state of the mo- 
tion \a = vTo)^!, with r]x = 0.1 (lower solid line), 0.15 (upper 
solid line). The dashed curve is given by \/T0exp(— Ft). Other 
parameters are k = 1, Ux = 5cA = 10, Qq/Aqa ~ 0.2, and 
flxgoSA/AoA = 0.1 (corresponding to F = 0.01). Figures (b) 
and (c) show the cavity field amplitude |(o(t))| (xlO) and the 
motional mode amplitude |(6a:(t))| for r]x = 0.1 and ijx = 0.15, 
respectively. Note that some of the structure (i.e., what looks 
like "beats") in |(d(t))| is due only to the finite number of 
points plotted. 

FIG. 5. Fidelity f{t) as a function of time for the pa- 
rameters of Figure 4 with ijx =0.1 (circles), and i]x = 0.15 
(squares) . 

FIG. 6. Cascaded atom-cavity systems for motional state 
transfer. The coupling between cavities is assumed to be uni- 
directional (facilitated, for example, by Faxaday isolators). 
During an ideal transfer no photons are detected by the pho- 
todetector, which monitors the output from the second cavity. 



FIG. 7. Schematic for state transfer between light and 
collective motional modes of a string of trapped ions, (a) 
First, a single- ion interaction with the cavity field mode is 
used to "receive" an incoming quantum state, after which (b) 
auxiliary lasers couple the single-ion and collective vibrational 
modes. 
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